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We study experimentally the motion of an intruder dragged into an amorphous monolayer of 
horizontally vibrated grains at high packing fractions. This motion exhibits two transitions. The 
first transition separates a continuous motion regime at comparatively low packing fractions and 
large dragging force from an intermittent motion one at high packing fraction and low dragging 
force. Associated to these different motions, we observe a transition from a linear rheology to a 
stiffer response. We thereby call "fluidisation" this first transition. A second transition is observed 
within the intermittent regime, when the intruder's motion is made of intermittent bursts separated 
by long waiting times. We observe a peak in the relative fluctuations of the intruder's displacements 
and a critical scaling of the burst amplitudes distributions. This transition occurs at the jamming 
point characterized in a previous study 1 and defined as the point where the static pressure 
{i.e. the pressure measured in the absence of vibration) vanishes. Investigating the motion of the 
surrounding grains, we show that below the fluidisation transition, there is a permanent wake of free 
volume behind the intruder. This transition is marked by the evolution of the reorganization patterns 
around the intruder, which evolve from compact aggregates in the flowing regime to long-range 
branched shapes in the intermittent regime, suggesting an increasing role of the stress fluctuations. 
Remarkably, the distributions of the kinetic energy of these reorganization patterns also exhibits a 
critical scaling at the jamming transition. 
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Jamming occurs, when a system develops a yield stress 
in a disordered state O |3l |4] and has been reported in a 
wide class of systems such as colloids [5 , foams [6j , emul- 
sions [Tj , granular materials [H [9] as well as in various 
model situations [3l [TOl [Til US]- One possible mecha- 
nism for such a change between fluid and solid like be- 
havior is that rearrangements of particles become pro- 
gressively slower while the stress relaxation time grows 
dramatically. The dynamics becomes spatially heteroge- 
neous and temporally intermittent, while the stress re- 
sponse appears more and more heterogeneous. A strin- 
gent manifestation of such inhomogeneities is the "stick- 
slip" response observed when the system is driven close 
to yielding and flows in rapid bursts. However the inter- 
play between density fluctuations and stress relaxation is 
still poorly understood Questions of interest are as fol- 
low: What is the nature of the re- arrangement events? 
How do these events depend on the external load and 
packing fraction? Also from a more fundamental view- 
point, whether the emergence of a yield stress coincides 
with dynamical arrest is still a matter of debate. 

Microrheology is a promising technique providing local 
probes of the dynamics in complex fluid [13 . Studying 
the motion of an intruder embedded in the material of 
interest, one is able to investigate the microscopic ori- 
gins of the complex-fluid behavior and in particular the 
link between microscopic mechanisms and macroscopic 
properties as given by conventional rheology. Applying 
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a force to the intruder, one explores the non-equilibrium 
and usually non-linear response, providing detailed in- 
sight into the structure-dynamics relationship. Previous 
drag experiments in colloids [5|, foams [14 , static [15] 
and shaken [16 granular media as well as simulations of 
structural glasses [10 were focusing on the velocity de- 
pendence of the drag force: proportionality is found for 
loose enough systems, reminding Stoke's law, while an 
increasing yield stress appears for denser packing. Stress 
fluctuations have been studied in details in [17 , and spa- 
tial reorganizations in [18 , however a clear picture filling 
the gap between spatial fluctuations and rheological ob- 
servations is still lacking. 

In the present paper, we investigate the motion of an 
intruder dragged with a constant force, within an amor- 
phous monolayer of horizontally vibrated grains, a sys- 
tem for which the jamming transition has been clearly 
identified and characterized in terms of the critical be- 
havior of the dynamics in a previous study [1 . At moder- 
ate packing fractions, and comparatively high force, the 
intruder moves rapidly as soon as the force is applied. 
Above some threshold value of the packing fraction which 
increases with the applied force, the intruder exhibits an 
intermittent creep motion with strong fluctuations rem- 
iniscent of a "crackling noise" signal. Simultaneously, 
the force-velocity relation evolves from a linear rheology 
to a stiffer response, thereby suggesting to call "fluidisa- 
tion" this first transition. A second transition is observed 
within the intermittent regime, when the intruder's mo- 
tion is made of intermittent bursts separated by long 
waiting times. This transition is signed by a peak in 
the relative fluctuations of the intruder's displacements 
and a critical scaling of the bursts amplitudes' distri- 
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butions. This transition occurs at the jamming point 
(j)j characterized in [T] and defined as the point where 
the static pressure {i.e. the pressure measured when the 
vibration is switched off) vanishes. In [1 , the authors 
demonstrate that dynamical heterogeneities become crit- 
ical at the transition. Here we investigate the motion of 
the grains surrounding the intruder. Below the fluidiza- 
tion transition a wake of free volume is observed behind 
the intruder and the fluidisation transition is marked by 
the evolution of the reorganization patterns around the 
intruder, going from compact aggregates in the flowing 
regime to long-range branched shapes in the intermittent 
regime, suggesting an increasing role of the stress fluctu- 
ations. The distributions of the kinetic energy of these 
reorganization patterns also exhibits a critical scaling at 
the jamming transition. 

This paper is an extended version of a recently pub- 
lished letter [19 , in which both transitions have been 
reported. The purpose of the present paper is to pro- 
vide a comprehensive study of the displacement fields 
surrounding the intruder and to take this opportunity to 
provide details on our analysis procedures as well as an 
extended discussion of our results. The paper organizes 
as follows : the experimental set-up and protocols are 
described in section HJ In section HIl we first introduce the 



Camera 



raw dynamical quantities and the phase diagram (II A), 



then we characterize the fluidisation (IIB) and jamming 



(IIC) transitions. The dynamics around the intruder is 
analyzed in section [ITU both by the relation between the 
average flow and the spatial fluctuations (III A) and by 



the evolution of the averaged free volume around the in- 
truder (IIIB). Finally, we study the bursts statistics in 



the intermittent regime close to jamming in section [TVj 
A general discussion and a few concluding remarks are 
given in section |V| 



I. EXPERIMENTAL SET-UP AND PROTOCOL 

The experimental set-up has been described else- 
where [Ij and we shall only recall here its most impor- 
tant characteristics and the modifications induced by the 
dragging procedure. The system is made of a mono- 
layer of 8500 bi-disperse brass cylinders of diameters 
dsmaii = 4 ± 0.01mm and di^ig = 5 ± 0.01mm laid out 
on a horizontal glass plate vibrated in its plane at a fre- 
quency of lOHz and with a peak-to-peak amplitude of 
10mm. The grains are confined in a cell, fixed in the lab- 
oratory frame, the volume of which can be adjusted by a 
lateral mobile wall controlled by a fim accuracy transla- 
tion platen, which allows us to vary the packing fraction 
(j) of the grains by tiny amounts ~ 5 x 10""^). The 

pressure exerted on this wall is measured by a force sen- 
sor (see fig{l]). The intruder consists in a larger particle 
) of same height introduced in the 
system and pulled by a mass via a pulley perpendicularly 
to the vibration. In all data presented here the resultant 
motion is strongly overdamped and the applied force can 
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FIG. 1: Top Experimental Setup (see text for details). Bot- 
tom Two samples of displacement fields around the intruder 
at for similar x-displacements of the intruder (Ax = 2.14 
and 2.08, in respectively 14 and 131 time steps). The arrows 
show the unmagnified grain displacements, while the interpo- 
lated amplitude field is in grey scale. 



be considered as constant. We use a fishing wire that 
stands over the other grains and doesn't disturb their 
dynamics. The time unit is set to one plate oscillation 
while the length unit is chosen to be the diameter of the 
small particles. The drag forces F are expressed as the 
ratio of the applied mass onto the total mass of grains in 
the ceh {Mtotai = 2.36bkg). 

We have been using two protocols, setting either the 
drag force or the packing fraction constant. We explore 
the force / packing fraction diagram along constant force 
lines using three different drag forces (Fi = 0.029, F2 = 
0.064 and F3 = 0.113) and varying the packing fractions 
from (j) = 0.8223 to 0.8418, and along constant packing 
fraction lines using three different packing fractions (^1 = 
0.8383, ^2 = 0.8304 and ^3 = 0.8399) and varying the 
drag force from F = 0.029 to 0.617, as sketched in fig|2] 

Starting from a low packing fraction ^, we gradually 
compress the system until it reaches a highly jammed 
state following the same protocol as in [1 . Then we 
stepwise decrease the packing fraction. In the absence 
of the intruder, it was shown on the one hand that the 
average relaxation time increases monotonically with the 
packing fraction and on the other hand that the dynam- 
ics exhibits strong dynamical heterogeneities, the length- 
scale and time-scale of which exhibit a sharp peak at an 
intermediate packing fraction. The authors have shown 
that the spatial correlations of these dynamical hetero- 
geneities exhibits a critical scaling at the transition. The 
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FIG. 2: Experimental conditions in a force - reduced packing 
fraction e diagram. Each point corresponds to a trajectory of 
the intruder into the media. The horizontal and vertical dot- 
ted lines correspond to the experimental exploration paths, 
either at constant force of packing fraction. O (resp. (8)) de- 
note the fluidized (resp. intermittent) behavior (see text for 
definition) . 



pressure measured at the wall in the absence of vibra- 
tion falls to zero precisely below that packing fraction, 
hence called the jamming transition (j)j. It is important 
to mention that for such high packing fractions, the struc- 
ture as given by the neighborhood relation among the 
grains is frozen on experimental time-scales. Hence the 
observed transition is to be understood as the jamming 
of a given frozen configuration. Accordingly the value of 
(j)j may slightly change from one run to another, since it 
may differ from one frozen configuration to another. In 
the present case, by monitoring the pressure at the wall 
while interrupting the vibration and without drag force, 
we could localize the jamming transition and obtain three 
close but different values : 0.8369 (Fi), 0.8383 (F2), 
0.8379(^3) and 0.8388 in the run where we explore iso- 
(j) lines. These values also are slightly smaller than the 
value = 0.8417 reported in pLj, maybe an effect of the 
geometrical distortion induced locally by the size of the 
intruder, twice larger than the other grains. In the fol- 
lowing we will use either the packing fraction (j) or the 
reduced packing fraction s = {(j) — <pj)/(j)j. 

In the present study, before each step, the intruder is 
removed from the position it has reached, replaced by 
one big and two small grains, and inserted at its initial 
position in place of one big and two small grains. Then 
a downward step in packing fraction is eventually done, 
and the system is kept under vibration up to one hour in 
order to "equilibrate" the configuration. At that point. 
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FIG. 3: Top left: Trajectories of the intruder for several 
packing fractions, from = 0.8306 (red) to = 0.8418 (blue), 
at constant force (F = F2). The trajectories at the highest 
packing fractions are truncated. Top right: log^Q of the av- 
erage displacement of the intruder Ax in function of the lag 
time r for the same packing fractions and the same drag force. 
Bottom left log^^Q of the average velocity of the intruder in 
function of the packing fraction 0, same packing fractions 
and drag force. Bottom right: Instantaneous velocities of 
the intruder at a high drag force (F = 0.363) for two very 
close packing fractions on both sides of the fluidisation and 
jamming transition, here indistinguishable. At = 0.8383 the 
intruder is always moving while at (f) — 0.8394 one can observe 
intermittent bursts of activity separated by long waiting in- 
tervals. 



the pressure has the same value as without the intruder 
for the corresponding packing fraction, indicating that 
the system has recovered from the small perturbation in- 
duced by the intruder's "teletransportation" . Only then 
the force is applied and the intruder is dragged through 
the cell, while its stroboscopic motion together with that 
of a set of 1800 surrounding grains in the center of the 
sample is tracked by a digital video camera triggered in 
phase with the oscillations of the plate. 



II. FLUIDISATION AND JAMMING 

A. Phase diagram 

When looking at the trajectories of the intruder along 
the drag direction x (see figjsj-top left), one immediately 
notices that the typical velocity dramatically changes 
within a tiny variation of the packing fraction : the 
intruder browses the entire system in a few time steps 
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FIG. 4: Left Ratio of negative displacements %6x~ in function of the reduced packing fraction e, for three different forces : 
Fi (♦), F2 (■) and F3 (a). For clarity, only one point with a ratio of is shown. Inset Fluidisation force in function of the 
reduced packing fraction e = (0 — Each point correspond to the moment where the ratio of negative displacements 

reaches 0. These points are also marked with arrows in the left and middle figures. Middle Ratio of negative displacements 
in function of the drag force F, for three different packing fractions : 01 ( ), 02 (o) and 03 ( ). Right Applied drag force in 
function of the average velocity of the intruder in a lin — lin plot at 0i. Inset Idem in a lin — log plot for different packing 
fractions : 0i ( ), 02 (o) and 03 ( ). Fits are eye-guiding affine (dotted) or logarithmic (plain) behaviors. 



for low packing fractions, and conversely seems to be 
arrested for the highest values of 0. For a given drag 
force, the average displacement probed over a lag time r, 
Ax(r) = (x(t + r) — x{t))t is roughly V.r, where V is the 
intruder's average velocity. V spans four decades from 
5 X 10~^ to 5 while varying the packing fraction of only 
a few percent (^0/0 = 2 x 10~^), illustrating the dra- 
matic freeze of the dynamics (figjsj-top right). For the 
highest packing fractions, one may notice a systematic 
bending of the curves at short times indicating that the 
intruder does not feel instantaneously the bias induced 
by the drag force. 



A closer inspection of the dynamics reveals two salient, 
distinct regimes of the intruder's motion. For loose pack- 
ings and large drags, the intruder moves continuously, 
while for dense packings and small drags the intruder's 
motion is highly intermittent. Anticipating on the follow- 
ing, let us call this transition "fluidisation" and empha- 
size that it is distinct from the jamming one as illustrated 
on the phase diagram (figj2|. Only for the largest forces, 
both transitions become asymptotically close, as illus- 
trated on figure [3}bottom left, where the instantaneous 
displacement of the intruder is shown for two packing 
fractions just below and above the transition. Just be- 
low the transition the intruder is moving continuously 
fast like in a fluid - note the time-scale on the horizontal 
axis- while just above it one observes violent bursts sep- 
arated by extremely long waiting times, indicating that 
the medium is not fluidized. We now characterize in more 
details the nature of these two transitions. 



B. Fluidisation 



The intruder's motion results from the drag force com- 
peting against the resistance of the surrounding grains. 
When the drag force is low enough - or if the packing 
is dense enough - the configuration can sustain the drag 
stress until some rearrangement of the force network, in- 
duced by the vibration, allows the intruder to move for- 
ward. In the meantime, the intruder's motion is cage- 
like and almost isotropic, going forward and backward 
roughly half of the time. Figure [4| left (resp. middle) dis- 
plays the percentage of time the intruder is going back- 
wards %5x~ as a function of the relative packing fraction 
e for the three constant forces Fi, F2, F3 (resp. as a func- 
tion of the dragging force for the three packing fractions 
01, 02, 03 )• In the extreme cases for which the system is 
most of the time stuck in jammed states %5x~ is very 
close to 0.5. As the drag force becomes stronger as com- 
pared to the resistance of the surrounding grains, the 
intruder will be less and less often blocked. As a result 
the percentage of backward steps will be smaller, and 
eventually will drop to when the intruder motion be- 
comes continuous. This is precisely how we have chosen 
to identify the fluidisation transition pointed out by an 
arrow on the figure. The most striking feature is that we 
could not observe fluidisation for the two packing frac- 
tions larger than 0j. %5x~ decreases with the force, but 
remains far from zero even for forces as large as 0.6, that 
is of the order of the force needed to drag all the grains 
on the glass plate in the absence of vibration. Our data 
suggests a divergence of the fluidisation line at the jam- 
ming transition : Ffiow ~ (see inset of figure [ij- left). 
Looking at the force-velocity relation as shown on fig- 
ure |4[right, we observe that in the intermittent motion 
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FIG. 5: Top left Pdf of the instantaneous displacements in 
the drag direction 6x at a constant force F — F2 for several 
values of cj). Inset same for the positive displacements only, 
in log — log scale. Top right Cumulated Pdf of the instan- 
taneous displacements in the drag direction 5x at a constant 
packing fraction cj) — for several values of F (right). Bot- 
tom Inverse average (left) and standard deviation over 
average cr^//i^ (right) of the instantaneous positive displace- 
ments in the drag direction as functions of the reduced pack- 
ing fraction e. Different curves correspond to different drag 
forces : Fi (green), F2 (blue) and F3 (red). 



regime, that is below the fluidisation line, F ~ hi{V) 
whereas F — Ffiow oc V in the fluidized phase. This 
demonstrates that the fluidisation line also separates two 
rheological behaviors. Such a fluidisation transition has 
been previously reported in other experimental studies. 
However a straightforward comparison can not be made 
without further enquiries, which we shall report to the 
discussion part. 
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Figure [sj-bottom left shows that the typical time 

the intruder takes to move one particle diameter, in- 
creases monotonically with s and faster than exponen- 
tially. The stronger the dragging force the sharper the 
increase. No signiflcant behavior is observed when cross- 
ing the jamming transition. On the contrary, flgure [5]- 
bottom right reveals a peak of cr+//i+ precisely at for 
the three dragging forces. Note that the peak sharpens 
when the drag force is stronger. Indeed, as discussed in 
the previous section, when F increases the fluidisation 
line approaches the jamming transition. This has two 
consequences. The displacement bursts become larger 
when they occur - as can be seen on flgure [5]-top left 
from the cumulated distribution of the intruder displace- 
ments for increasing forces - and the range of packing 
fraction separating the non fluctuating continuous mo- 
tion in the fluidized regime from the strongly intermit- 
tent one at jamming shrinks. Finally let us mention that 
in the intermittent regime the distributions of 6x~^ decay 
as power-laws with an exponent close to —3 (see inset 
of flgure [sj-top left), a result similar to the one reported 
in a recent simulation, where a probe is dragged into an 
assembly of harmonically repulsing disks [l^l . 

The above results reveal that the strong spatio- 
temporal heterogeneities of the dynamics reported in the 
absence of intruder [1 can also be seen in the response 
of the intruder to the dragging force. This is related to 
the fact that microrheology gives access to non-linear re- 
sponses and thereby probes the strength required to pull 
free the probe from the transient local structure. Before 
addressing in further details the critical nature of the 
fluctuations reported here above, we will concentrate on 
the response of the grains surrounding the intruder. This 
will ultimately allow us to perform the statistical analysis 
of the rearrangement events as a whole, involving both 
the intruder and the surrounding grains. 



III. AROUND THE INTRUDER 



C. Jamming 



A. Average flow and spatial fluctuations 



We now characterize the intruder dynamics when ap- 
proaching the jamming transition. When looking at the 
intruder's instantaneous displacements Sx along the drag 
direction in the intermittent regime (flgjsj-bottom right 
and flgj9}left), one immediately notices very strong fluc- 
tuations, with bursts of widely fluctuating magnitude in 
the direction of the drag. As a result the distributions of 
Sx (flgjsj-top left) exhibit an important skewness toward 
the positive displacements for the packing fractions above 
fluidisation. In order to characterize these positive dis- 
placements (5^, we compute the average value over time 



As already suggested in flgure [T] the instantaneous dis- 
placement fleld around the intruder during a burst is 
rather complex. They are typically asymmetric, with 
a main vortex on one or the other side of the intruder. 
Averaged over many bursts (see flgure [6]-top right), the 
displacement fleld recovers the intruder's left-right sym- 
metry and exhibit a two vortices pattern. Further in- 
sight into the dependence of this flow on the packing 
fraction is obtained when looking at the proflles of the x- 
component of the velocity along the direction of the drag 
(x-direction) and perpendicular to the drag (?/- direct ion) 
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FIG. 6: Left Velocity profiles : average velocity in the x- 
direction along the horizontal (top) and vertical (bottom) 
axis, for 15 packing fractions. Inset Same in absolute val- 
ues and log — lin. The dark dotted line is an exponential 
decaying guide with a characteristic length = 7. Top- 
right Displacement field around the intruder at = 0.8386, 
averaged over time. The arrows are magnified by a factor 3. 
Bottom-right Average percentage of particles in the con- 
nex clusters of particles faster than Vmaxf^ in function of the 
reduced packing fraction. Error bars show the standard devi- 
ation. The grey area correspond to the intermittent regime. 
In ah plots, F = F2. 



(see figure |6]-right top and bottom). As evidenced in the 
insets, one observes that the velocity field decreases ex- 
ponentially with the distance to the intruder, and amaz- 
ingly that the typical length-scales A^^ :^ 7 and A^ :^ 10 
associated with such a dependence are totally indepen- 
dent from the packing fraction. This is also confirmed 
by the localization of the center of the vortices, easily 
located on the ^/-profile, which remains at the same dis- 
tance from the intruder at all packing fractions. Hence, 
as far as the averaged flow is concerned, only its overall 
magnitude depends on the packing fraction and scales 
like y, the averaged velocity of the intruder. 

In order to study the fluctuations of the displacement 
field, one could imagine to subtract the averaged field; 
however without a prescription for the analytical form of 
this field or better statistics, this is a rather uncontrolled 
operation. An alternative procedure is to define first a 
maximal displacement for the grains, Vmax^ as the aver- 
age displacement of the intruder's closest neighbors, and 
then to define the "fast particles" cluster as the connex 
cluster of particles around the intruder which move more 
than Vmax/"^' The average number of particles in this fast 
cluster is plotted on figure |6]-bottom right as a function 
of the packing fraction: it grows from 10% to 80% in the 



tiny interval between (j) = 0.8306 and = 0.8409, under- 
lining how sharply the spatial extent of the reorganiza- 
tions grows as the system goes toward dynamical arrest. 
It is however surprising that the number of fast parti- 
cles increases, considering that the averaged flow scales 
entirely with the velocity of the intruder which we have 
scaled out by defining the "fast" particles relatively to 
the average velocity around the intruder. Such a differ- 
ence between the averaged flow behavior and the instan- 
taneous one can only be explained by the existence of 
strong heterogeneities in the instantaneous fields. The 
amplitude of the fluctuations of the number of fast par- 
ticles, as indicated by the error bars on figure |6]-bottom 
right is already an indication that it is indeed the case. 

A further characterization of these heterogeneities is 
provided by the shape of the clusters made of the p% 
fastest particles at each time step; typical examples with 
p = 15% are shown on figj7[top for three packing frac- 
tions below, at, and above ^0 the fluidisation packing 
fraction. One clearly observes that these instantaneous 
rearrangements evolve with the packing fraction from 
dense to branched patterns with long chains spanning 
the whole system. The number of neighbors rineigh inside 
these clusters is a good indicator of the level of branching 
since its typical value is 6 for perfectly dense clusters and 
2 for perfect strings. The well defined plateau in the de- 
pendence of (rineigh) with p (figure [Tj-middle center) cor- 
responds to the clusters which contain enough particles to 
have a non-trivial shape but do not reach the boundaries 
of the acquisition field. The average number of neighbors 
at the plateau (figj7[middle right) decreases significantly 
with the packing fraction providing a quantitative evalu- 
ation of the clusters' evolution from dense aggregates to 
branched structures. The contour fractal dimension (ij 
calculated with the classical compass method [21] also 
gives insight on branching. The mean value ((ij)(p) has 
a nice maximum for the range of p corresponding to the 
plateau in (rineigh) {p)- For a value of p = 40%, that 
is for clusters large enough to assess the fractal dimen- 
sion without spanning the acquisition field, we observe a 
significant increase of (ij((/)) at (/)o, the packing fraction 
where the intruder's motion becomes intermittent. 

Altogether, the patterns observed in the instantaneous 
displacement field during the relaxation events suggest 
that the force chains network, whose importance has been 
visually exemplified in photo-elastic disks experiments of 
an intruder dragged in a Couette cell [9 , starts playing 
a significant role as soon as the system crosses the flu- 
idisation transition. However, one can not elude the role 
of the density field, which we now investigate computing 
the free volume field around the intruder. 



B. Free volume 

We then examine the density field around the intruder. 
For that purpose, we compute the free volume field ex- 
tracted from Laguerre's tesselation [4? of the configura- 
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FIG. 7: Top Typical contours of clusters of the p = 15% 
fastest particles, for three packing fractions, from left to right : 
= 0.8337 (0 < 00), = 0.8358 (0 = 0o) and = 0.8396 
(0 > 00 ). The intruder is marked by a black disk and moves 
from left to right. Middle left Average number of neighbors 
{Tineigh) in fuuctiou of p, for different packing fractions. The 
black dotted line is at p = 40%, in the middle of the plateau. 
Middle right (rineigh) in the clusters of p = 40% fastest 
particles, as a function of e. The grey area correspond to 
the intermittent regime, and error bars show the standard 
deviation. Bottom left Average fractal dimension of the 
clusters contours (dj) as a function of p for different packing 
fractions. The black dotted line is at p = 40%. Bottom 
right (dj) as a function of £ for p = 40%. The grey area 
correspond to the intermittent regime, and error bars show 
the standard deviation. F = F2 for all plots. 



tions at each time step: the free volume V/ of a Voronoi' 
cell is the area difference with the minimal possible reg- 
ular hexagon around the grain normalized by the surface 
of the grain. Figure [8] displays the averaged free volume 
field around the intruder for three packing fractions, be- 
low, above, and at the jamming transition. Above 0j, 
there is a very small amount of decompaction around 
the intruder without significant signature of the intruder 
motion, seemingly an effect of the size difference between 
the intruder and the surrounding grains disturbing the 
local organization of the packing. On the contrary, as 
the packing fraction decreases below jamming, two holes 
grow, first on the back of the intruder and then on its 
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FIG. 8: Top and bottom- left Average free volume fields 
around the intruder for three packing fractions: = 0.8327 
(0 < 0j, top- left), = 0.8386 (0 ^ 0j, top-right), = 0.8405 
{(f) > (f)j, bottom-left). Bottom-right Average free volume 
(V/) before (v) and after (a) the intruder in function of the 
reduced packing fraction. Error bars represent the standard 
deviation over time of the free volume in each point of space, 
averaged on the computation window. 



front. 

As a consequence, one observes a clear signature of the 
intruder's motion, namely the apparition of a forward- 
backward asymmetry which can be quantified by com- 
puting for instance the average free volume in front of 
and behind the intruder (figure [sjbottom right) [43^. One 
observes that above , there is no more sign of the asym- 
metry associated with the intruder motion. 

Altogether, we have seen that the averaged flow around 
the intruder is surprisingly insensitive to the transition, 
apart from a simple scaling factor indexed on the average 
velocity of the probe. On the contrary, the fluctuations 
reveal a rather complex interplay of the density and stress 
fluctuations. Below the fluidisation transition, the in- 
truder motion is dominated by important free volume re- 
arrangements, which concentrate in dense clusters of fast 
moving grains around the intruder. As soon as the sys- 
tem enters the intermittent regime, these clusters start to 
branch, indicating the existence of inhomogeneities in the 
rigidity of the system. The instantaneous displacements 
spread on larger and larger scales; and progressively the 
free volume fluctuations become distributed throughout 
the system. Above these last type of rearrangement 
completely dominate the intruder motion. 
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IV. BURST STATISTICS CLOSE TO JAMMING 

In the above section, we have seen that the tempo- 
rahy intermittent motion of the intruder also corresponds 
to spatiahy heterogeneous rearranging regions. We wih 
now investigate the intruder displacements and the ki- 
netic energy of these regions, focusing on the tempo- 
ral correlations at the root of this dynamics of bursts. 
Figure [9|bottom displays the instantaneous displacement 
5x{t), and the kinetic energy involved in a rearrange- 
ment dE{t) = H^xf (t), where the sum is performed on 
the clusters of particles faster than Vmaxf^^ as defined in 



section III A For three different packing fractions cho- 
sen above the fiuidisation transition, one clearly observes 
a sequence of distinguishable and well separated bursts, 
suggesting "crackling noise" -like signals [22, 23 . 

Note that the signals of 6x{t) and SE{t) are very sim- 
ilar. However the correlation is not as strong as one 
would first imagine by visual inspection. Each burst in 
the displacements signal corresponds indeed to a burst 
in the kinetic energy one but the bursts might be shifted 
in time and the intruder's displacements sometimes pre- 
cede, sometimes follow the energy ones. In addition, the 
amplitude of the bursts only weakly correlate and there 
are some occasions where there is a burst of activity in 
the system without a displacement of the intruder itself 
(see an example on fig.[9|, but these are marginal events. 

The first step in analyzing such kind of signals is to 
coarsegrain the statistical distribution of the jumps to 
capture the size and duration of the relevant bursts. One 
common way to do so is to define a threshold, or reso- 
lution coefficient, in order to delimit the temporal limits 
of each burst : given a signal u{t) and a threshold one 
can define bursts or avalanches as the connex portions of 
the signal that stand above u ; every burst have a du- 
ration and an integrated amplitude S^^ (see figjoj-top 
left). Choosing the threshold is a delicate problem given 
the large variations of the displacements' amplitudes, the 
complete loss of small events for the loosest packings and 
the fact that the background noise limits the detection of 
small bursts at high packing fractions. In the following 
we have set all the thresholds at the value where the dif- 
ference between the cumulated distributions of the local 
minimums and maximums of u{t) exhibits a peak. On 
one hand this corresponds to the point where the mini- 
mums and maximums in u{t) are best separated, and on 
the other hand it points out to the threshold value for 
which the number of bursts is maximal, naturally enforc- 
ing the statistics. We have checked that multiplying the 
so-obtained thresholds by a factor from 1/2 to 2 doesn't 
change the following results. 



First we observe on figure [10}top right that the waiting 
times separating two succesive bursts are exponentially 
distributed the characteristic time of which r^^ :^ 2 is in- 
dependent of the packing fraction (see inset). Two con- 
clusions can be drawn: (i) the bursts follow a Poissonian 
process and can be considered as independent events, and 
(ii) Tw remains very small compared to the total time of 
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FIG. 9: Top left Definition of the bursts of activity : given 
a threshold on a signal each connex portion of signal above 
the threshold defines a burst. Tu is the duration of the burst 
while Su is its area (in grey). The waiting times between 
bursts are noted ■ Top Right Cumulated probability den- 
sity functions of the time intervals between bursts of dis- 
placements for 15 packing fractions from = 0.8306 (red) 
to 0.8418 (blue), at the constant drag force F2. Inset Aver- 
age waiting time Tw for the same packing fractions estimated 
by an exponential fit. Error bars represent the 95% confi- 
dence bounds of the fit. Bottom left Instantaneous displace- 
ments Sx{t) of the probe for three packing fractions. Bottom 
right Energy of the grains 5E{t) around the intruder for the 
same packing fractions. Note that all peaks can be found in 
both signals, except some rare events like the peak number 4. 
F — F2 for all bottom plots. 



the experiments. We catch a large number of bursts, 
even at the highest packing fractions. Note that the use 
of a constant threshold to define the bursts would have 
lead to a dramatic increase of these waiting times in the 
intermittent regime and for the highest packing fractions 
the number of bursts would have been vanishingly small. 

We focus now on the statistics of the bursts themselves. 
Let us note T and L, respectively Te and the duration 
and the integrated amplitude of the bursts recorded on 
the signal of the displacements of the probe, respectively 
of the kinetic energy of the surrounding fastest grains. 
Obviously these quantities are not independent one from 
another. The first step in the analysis (see figure [lO]- left), 
shows that : 



L(T,(/)) =Lo((/>)Ti/^ and E{Te 



Eom'J\ (2) 



where Lq and £^0 can be interpreted as the typical dis- 
placement and energy associated with a burst and z = 
2/3 is usually called the dynamical exponent. Lq depends 
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FIG. 10: Top left: Rescaled length L/Lq as a function of T. Inset Scaling factor Lq as a function of e. The black dotted 
line corresponds to the average speed of the intruder. Top middle: Cumulated Pdf of the burst durations T in the intruder's 
displacement signal. Inset Cut-off 0(e) of the distribution of T. The eye-guiding lines are power laws with an exponent —2/3. 
Top right: Cumulated Pdf of the rescaled length L/Lq. Inset Cut-off ^{e) of the distribution of L/Lq. The eye-guiding lines 
are power laws, with an exponent —1. Bottom left: Rescaled energy E/Eq as a function of Te- Inset Scaling factor Eq 
as a function of e. Bottom middle: Cumulated Pdf of the burst durations Te in the energy signal. Inset Cut-off Oe{£) of 
the distribution of Te- The eye-guiding lines are power laws with an exponent —2/3. Bottom right: Cumulated Pdf of the 
rescaled energy E/Eq. Inset Cut-off T(£) of the distribution of E/Eq. The eye-guiding lines are power laws, with an exponent 
-1. 



on ^ in a similar way as the average velocity of the in- 
truder does, consistently with the previous observation 
that the average waiting time is independent from cj). 



relations: 



The second observation concerns the distributions of 
T, Te and L/ Lq^E / Eq. Again we find identical behaviors 
for the intruder displacement and the kinetic energy of 
the fastest surrounding grains (fig. [Toj- middle and right). 
All the above quantities are largely distributed, with a 
large value cutoff, which depends on the packing fraction. 
Visual inspection of the tails of the distribution already 
indicates that the cutoff dependence on the packing frac- 
tion is not monotonic. In all plots the jamming packing 
fraction corresponds to the light green curves. The distri- 
butions are the largest for that precise packing fraction. 
Such behaviors can be encoded in the following scaling 



p{T) (X T-^i+")./ 



p(Te) a T-^'^-\fE 
p{L/Lo) a (L/Lo)-(^+^) 



p{E/Eo) a {E/Eo) 



Te 

L 

E 

9e 



(3) 
(4) 

(5) 
(6) 



where the exponents a 1/2 and :^ 1/3 satisfy the 
expected relation a.z = jS. Note that given the lack of 
statistics these exponents are not the results of a fit but 
only indicative values. One then easily extracts the scal- 
ing variables, which are the cutoffs of the distributions 
and are simply proportional to their mean values [44]. 
Their dependence on the packing fraction are plotted in 



the insets of figure 10 and exhibit a sharp peak at the 
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FIG. 11: Comparison of the evolution of dynamical length- 
scales : the cut-off of the distribution of the bursts of displace- 
ment's amplitudes ^ is plotted for 3 different forces (O Fi, □ 
F2 and A F3) and compared to (•) the long-range correlation 
length ^4, revealed from the dynamics of the grains without 
the intruder. The dotted and plain curves are guiding the eye 
like 



jamming transition, suggesting a critical behavior and 
thereby ensuring the self-consistency of the above scal- 
ing analysis. This critical behavior is well described by: 



(9(6:) (x6>£;(£:) (xe-^ 
i{e) (X T(£:) (X e'^ , 



(7) 
(8) 



where r] 
T] = u.z. 



-2/3 and :^ — 1 also satisfy the relation 



Altogether the above analysis provides strong evi- 
dences of a critical behavior of the intruder motion at 
the jamming transition, enforcing the first evidences of 
dynamical criticality reported in 
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reports 



[1 . Figure 

the lengthscale ^ measured here for the three different 
dragging forces, together with the dynamical correlation 
length ^4, which is the correlation length of the local 
density relaxations, and which was computed in pT. The 
authors reported a dependence of ^4 with the distance to 
jamming compatible with ^4 ~ however it can not 

really be discriminated from ^4 ~ e"^, as reported in the 
present work for the lengthscale ^. The fact that both 
lengthscales behave in a similar way at the transition sug- 
gests a fluctuation-dissipation-like relation between the 
non-linear response studied here and dynamical hetero- 
geneities, as recently discussed in the context of Mode 
Coupling Theory [24 . In the following we discuss our 
results as compared to other microrheological studies in 
dense systems of particles, and try to precise the kind of 
criticality we are facing at jamming. 



V. DISCUSSION AND CONCLUDING 
REMARKS 

We have experimentally studied the motion of an in- 
truder dragged into an amorphous monolayer of horizon- 
tally vibrated grains at high packing fractions. A first 
"fluidisation" transition separates a continuous motion 
regime, where the force- velocity relation is afiine, from an 
intermittent motion one, where the force depends loga- 
rithmically on the velocity. The force threshold increases 
with the packing fraction and seemingly diverges at the 
jamming transition defined as the packing fraction where 
the pressure goes to zero in the absence of vibration. Be- 
low this threshold the intruder motion is intermittent. 
We have reported the existence of a second transition, 
where the fluctuations associated with this intermittent 
motion exhibit a critical behavior. This evidence is sup- 
ported by the analysis of the motion of the probe but 
also by that of the surrounding grains. The bursting 
events are characterized by increasingly heterogeneous 
patterns in the instantaneous displacement field around 
the intruder and an increasingly spreading redistribution 
of the free volume. 

Several other experimental, numerical and theoretical 
investigations of the response of a locally driven particle 
in a dense system report similar observations. Our aim 
here is to discuss similarities and differences among these 
studies, in an attempt to shed some light on sometime 
apparently controversial results. To our knowledge the 
first measurements were conducted in colloidal suspen- 
sions [5 . At low packing fractions, the authors report a 
linear force- velocity relation. At higher packing fraction, 
but below the glass transition a yield stress Fq devel- 
ops below which the probe remains trapped. For larger 
forces, the probe is delocalized by the applied force, the 
bead is pulled with a fluctuating velocity and a non-linear 
force- velocity relation holds. As emphasized by the au- 
thors themselves, the existence of a force threshold below 
the glass transition, where the spontaneous fluctuations 
alone still allow the particles to escapes their cages and 
relax, is rather intriguing. Indeed, in a recent work [25] , 
the force-velocity relation in dense suspensions has been 
investigated theoretically in the context of Mode Cou- 
pling Theory and compared to some of the above ex- 
perimental data p and to numerical simulations of a 
slightly polydisperse quasi-hard-sphere system undergo- 
ing strongly damped Newtonian dynamics. This time a 
force threshold is predicted to delocalize the probe parti- 
cle above the glass transition. Below the glass transition 
a strongly non-linear force-velocity relation is predicted 
and reported in the simulations. This steep dependence 
of the velocity on the applied force could explain the 
above observation of a threshold in an experiment where 
the lower resolution on the velocity is always finite. 

Recent simulations consider the motion of a single 
probe particle driven with a constant force in a binary 
mixture of two-dimensional disks with stiff spring repul- 
sive interactions at zero temperature. |20l |26]. As the 
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packing fraction is increased towards jamming, the av- 
erage velocity of the probe particle decreases and the 
velocity fluctuations show an increasingly intermittent 
or avalanche-like behavior. The velocity distributions 
are exponential away from jamming and have a power 
law character when approaching it within less than one 
percent. These observations are very similar to those 
reported in the present study. The averaged velocity 
dependence on the packing fraction exhibits the same 
faster than exponential decrease when approaching jam- 
ming and the power-law character of the velocity distri- 
butions close to jamming are identical. However the au- 
thors report the existence of a critical threshold force that 
must be applied for the probe particle to move through 
the sample, which increases when increasing the pack- 
ing fraction a 60 i;e jamming, whereas we never observed a 
complete arrest of our probe. In the simulations, once the 
probe stops, it cannot move any more because there is 
no temperature. In the present experiments the vibration 
allows the system to explore new force chains configura- 
tions and thereby to provide the intruder new opportu- 
nities of moving. As a final remark on the force-velocity 
relation, we would like to emphasize that the affine rela- 
tion, which we observe above threshold when the motion 
is continuous, is not to be confused with the linear re- 
sponse of Stoke's law. On the contrary, we believe that 
our observations correspond to the highest force regime 
reported in [25] , where the non-linear response ultimately 
recovers a linear behavior. 

An important difference between our system and the 
systems described above is that we are dealing with fric- 
tional particles. The role of friction close to jamming 
is an important issue. Recent important progress have 
been made in this matter [271 EH El], but many open 
questions remain unsolved. The isostatic criteria for me- 
chanical stability z = 2d valid for frictionless systems a 
priori turns into a double inequality d-\-l < z < 2d in the 
presence of friction, suggesting that friction could blur 
the critical character of jamming reported in frictionless 
systems. It was suggested in [20 that the absence of 
power-law behavior in the velocity distribution observed 
in [17] might indeed be an effect of friction. Our ob- 
servations clearly demonstrate on the contrary that the 
jamming transition remains critical in the presence of 
friction. Such a result suggests the existence of a general- 
ized isostaticity criteria for frictional systems as proposed 
recently [29 . 

Not only the intruder motion exhibits the strongest 
intermittency at jamming, but also the statistics of the 
bursts obey critical scalings. The later are reminis- 
cent of many other phenomena such as earthquakes [30] , 
Barkhausen noise [23 , crack tips dynamics in heteroge- 
neous materials [31 , generically assimilated to the so- 
called crackling noise [22^. It is also a distinct behavior 
of various random failure and load redistribution mod- 
els [32], which can be used to describe stress redistribu- 
tion in stick-slip granular experiments [33]. These simi- 
larities, even at some quantitative level - see in table [ij 



System [Ref.] 


1 + a 


1 + /3 


z 


Here 

Barkhausen noise |23j 
Crack tips ^35^, ^36^, '37', '38; 
Stick skip (grains) 39, 40, 41 


- 3/2 
1.5 - 2.3 
1.2 - 1.9 
1.8 - 2.1 


^1/3 
1.24 - 1.77 
1.47 - 1.51 
1.67 - 1.94 


- 2/3 
1/2 - 2/3 

/ 
/ 



TABLE I: Comparison of the critical exponents for several 
systems displaying "crackling noise" -like signals. 



a comparison of some of the exponents obtained in dif- 
ferent systems - should not hide important differences 
among these phenomena even at the conceptual level. 
First the kind of criticality reported here has little to 
do with self organized criticality, often associated with 
crackling noise observations [22 . Second in many situa- 
tions, an external parameter is increased (i.e. the loading 
force, the magnetic field) and the system fails once it has 
overcome some randomly distributed threshold. A con- 
trario in the present situation, the dragging force is kept 
constant, the system is vibrated and it explores succes- 
sive configurations. To what extent the study of these 
different dynamics are complementary is an interesting 
issue. Further investigations with an intruder dragged at 
constant velocity should provide interesting clues in this 
direction. 

Altogether, investigating the motion of a probe 
dragged at constant force in a dense granular media, we 
have identified a force threshold diverging at the jam- 
ming transition, below which the motion of the probe is 
intermittent and exhibits criticality at jamming. More 
generally, we believe that micro and macro-rheological 
studies combined to statistical observations such as dy- 
namical correlations are key elements to further investi- 
gate the underlying mechanisms of the jamming transi- 
tion in frictional systems. 
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